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Recent developments on the determination of the spin- 1/2 spectrum of the nucleon in full QCD 
are presented. Our focus is on the PACS-CS 2+1 flavor configurations made available through the 
ILDG. Using correlation matrix techniques, in which a wide variety of gauge-invariant Gaussian- 
smeared fermion-propagator sources and sinks are considered, excited states are determined. We 
consider several correlation matrices of various sizes, each constructed with a different set of 
basis interpolators, in order to demonstrate the invariance of the eigenstates on the basis choice. 
Of particular interest is the approach to the elusive Roper resonance and we report preliminary 
results in full QCD. 
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1. Introduction 

The first positive-parity excited state of the nucleon, known as the Roper resonance, N\ (1440 
MeV) Pn, has presented a long-standing puzzle since its discovery in the 1960's due to its lower 
mass compared to the adjacent negative parity, Nj (1535 MeV) Sn, state. In constituent quark 
models with harmonic oscillator potentials, the lowest-lying odd-parity state naturally occurs below 
the Pn state [jl|, Qj. In nature the Roper resonance is almost 100 MeV below the Sn state. 

Lattice QCD is very successful in computing many properties of hadrons from first principles. 
In particular, the ground state of the hadron spectrum is a well understood problem. However, the 
excited states still provide a significant challenge. The first analysis of the positive parity excitation 
of the nucleon was performed in Ref. []|] using Wilson fermions and an operator product expansion 
spectral ansatz. Since then several attempts have been made to address these issues in the lattice 
framework. 

The "variational method" [Q, ||] is one of the state-of-the-art approaches to hadron spec- 
troscopy, which is based on a correlation matrix analysis. The identification of the Roper state 
with this method wasn't successful until recently. In Refs. ^ [7|] a low-lying Roper state has been 
identified with this approach by employing a diverse range of smeared-smeared correlation func- 
tions. Also, with this method, a physical level ordering between the Roper and ground state 
is observed in Ref. 

Recent developments of algorithms and computational power have enabled us to explore this 



physics in full QCD. Some recent full QCD results has been presented in Refs. H9L MOL 11]. In 



this paper, we present preliminary results for the excited states of the nucleon using the PACS-CS 



2+1 flavor configurations [12]. The results are presented from the variational method using various 



fermion source and sink smearings to construct correlation matrices. 

2. Variational Method 

The two point correlation function matrix for a p = baryon can be written as 

Gfj(t) =£Tr sp {r ± {Q|^(x)^(0)|n)} 5 (2.1) 

X 

= Y^K^ e ~ mat , (2-2) 
a 

where Dirac indices are implicit. Here, A" and X" are the couplings of interpolators %i an d Xj at 
the sink and source respectively and a enumerates the energy eigenstates with mass m a . T± = 
5(70 ± 1) projects the parity of the eigenstates. 

Since the only t dependence comes from the exponential term, one can seek a linear superposition 
of interpolators, Xj u p sucn that, 

Gijfo + AO uj = e- m " At Gij(to) uj , (2.3) 



for sufficiently large to and to + At. More detail can be found in Refs. [13, 14, 15]. Multiplying the 



above equation by [Gjj(to)] 1 from the left leads to an eigenvalue equation, 



[(G(to))- l G(t + At)] ijU f =c a uf, (2.4) 
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where c c 



-m a At 



is the eigenvalue. Similar to Eq. (2.4), one can also solve the left eigenvalue 



equation to recover the v a eigenvector, 



vf[G(t + At)(G(t Q )) 



-l 



■■-r a v a 



(2.5) 



The vectors uj and vf diagonalize the correlation matrix at time to and to + At making the projected 
correlation matrix, 



The parity projected, eigenstate projected correlator, 

is then analyzed using standard techniques to obtain the masses of the states, a. 



(2.6) 



(2.7) 



3. Simulation Details 



The PACS-CS 2 + 1 flavor dynamical-fermion configurations [ 12 ] are used. The non-perturbatively 



& id) -improved Wilson fermion action and Iwasaki-gauge action [ ]16[ ] are employed. The lattice 
volume is 32 3 x 64, with j8 = 1.90 providing lattice spacing a = 0.0907 fm and csw = 1-715 [JO]. 
Five kappa values with degenerate up down quarks have been considered, i.e. k uc / = 0.13700, 
0.13727, 0.13754, 0.13770, 0.13781, with k s = 0.13640. In this paper, we also consider the Som- 
mer scale [18], with a as described in Ref. [12]. The results are presented for k ui j = 0.13754,0.13770, 
for which ensembles of 350 configurations are considered. As with the quenched case & various 
levels of gauge invariant Gaussian smearing JT9|] are applied at the fermion source (t = 16) and at 



the sink. We consider 4, 9, 16, 25, 35, 50, 70, 100, 125, 200, 400, 800, 1600 sweeps, corresponding 
to rms radii in lattice units of 1.20, 1.79, 2.37,2.96, 3.50,4.19,4.95,5.920, 6.63,8.55, 12.67, 15.47, 
16.00. The error analysis is performed using the jackknife method, where the # 2 /dof is obtained 
via a covariance matrix analysis. Our fitting method is discussed extensively in Refs. [^, |l4| ]. 
The nucleon interpolator we consider here is the local scalar-diquark interpolator [f|, [20|], 



-abc • I ,Ta 



u la {x)Cy 5 d h {x))u c (x). 



(3.1) 



4. Results 

We consider several 4x4 correlation matrices. Each matrix is constructed with different 
sets of correlation functions, each set element corresponding to a different number of sweeps of 
gauge invariant Gaussian smearing at the source and sink of the correlators. This provides 
a large basis of operators as described in Table [T], providing a wide range of overlap with energy 
eigenstates. 

Consistency of the extracted masses is manifest in Fig. [j]. In particular, the ground and Roper 
states are robust. The highest energy state shows some basis dependency (smearing dependency) 
which is to be expected as this state must accommodate all remaining spectral strength 
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Table 1:4x4 correlation matrix bases. 



Sweeps — > 


16 


25 


35 


50 


70 


100 


125 


200 


400 


800 


Basis No. | 


Bases 


1 


16 




35 




70 


100 
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16 




35 




70 




125 
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16 




35 






100 




200 






4 


16 




35 






100 






400 




5 


16 






50 




100 


125 








6 


16 
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Figure 1: (Color online). Masses of the energy states for various 4x4 correlation matrices (bases), 
given in Table [j], for k ul \ = 0. 13770, over 50 configurations. 



It is noted that basis operators that are linearly dependent will cause the eigenvalue analysis 
to fail as there will be a singularity in the correlation matrix. The fact that our analysis succeeds 
indicates that our choice of operators access an equal number of dimensions in the Hilbert space. It 
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is interesting to examine the stability of the masses to different choices of bases to ascertain whether 
one has reliably isolated single eigenstates of QCD. The relevant issues are: (i) whether or not the 
operators are sufficiently far from collinear that numerical errors do not prevent diagonalisation of 
the correlation matrix and, (ii) whether or not the states of interest have significant overlap with 
the subspace spanned by our chosen sets of operators. Since our correlation matrix diagonalisation 
succeeded, except at large Euclidean times where statistical errors dominate, we conclude that our 
operators are sufficiently far from collinear. 

Basis numbers 4, 7, 9 and 10 contain higher smearing-sweep counts of 400 and 800, which 
results in a significant enhancement of errors for the second and third excited states. It is noted that 
the sources with sweep counts of 400, 800 and 1600 are very challenging as the smearing radii for 
these sources are close to the wall source. The poor signal-to-noise ratio for these sources make 
the correlation matrix analysis more challenging and the eigenvalue analysis becomes unsuccessful 
for a large number of variational parameters (to, At). Therefore, the sources 400, 800 and 1600 are 
undesirable to work with. 



> 
O 

S2 







1 1 1 1 


I 1 1 1 1 


1 1 1 1 1 


1 , , , 


' 1 

O 


1 1 1 

g.s 










A 


1 st e.s : 










□ 


2 nd e.s" 






3 5 I 




O 


3 rd e.s - 


15 £ 5 5 






si 






; s a a 2 


ssss 


a a a 








-0000 


0000 


000 


O O 






■ 1 1 1 


1 . , , , 


1 . , , , 


1 , , , 


■ 1 





17 



18 



19 



20 



21 



22 



Figure 2: (Color online). Masses of the nucleon, Nk + states, from the projected correlation functions as 



shown in Eq. (2.7). Each set of ground (g.s) and excited (e.s) states masses correspond to the diagonalization 
of the correlation matrix for each set of variational parameters to (shown in major tick marks) and At (shown 
in minor tick marks). Figure corresponds to k uc [ = 0.13754 and for the 3rd basis. 



The agreement among the three lowest lying eigenstates is remarkable and verifies that our 
approach successfully isolates true eigenstates Jo], [l4|]. 

Basis number 3 has good diversity including both lower and higher smearings which is neces- 
sary for the extraction of masses over the entire heavy to light quark mass range. As a result, basis 
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Figure 3: (Color online). As in Fig. ||, but masses are calculated from the eigenvalues as shown in Eqs. (2.4) 
and (O). 



number 3 is considered for the following more extensive correlation-matrix analysis. 

Consistency of the extracted masses from projected correlation functions over the variational 
parameters to and At is evident in Fig. ||, whereas a significant dependence of the extracted masses 
from the eigenvalues is evident in Fig. [| Asymptotically, these energies agree with those of Fig. ||. 
These aspects of the extracted masses on variational parameters is in complete accord with our 



quenched analysis presented in Refs. [pLM,^, 14] 



5. Conclusions 



In these proceedings, we have presented a systematic correlation matrix analysis method using 
a variety of fermion field smearings at the source and sink to extract excited-state energies of the 
nucleon. Of particular interest is the Roper resonance from 2+1 flavor QCD. The negative parity 
channel will also be investigated to obtain information about the level orderings between the Roper 
and Nl ground state in full QCD. Performing this analysis technique at all the quark masses will 
be the subject of future investigations. 
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